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Longitudinal  spatial  effects  in  the  nonlinear  medium  are  analyzed  in  a  model  for  mirrorless  optical  bistability  due  to 
induced  absorption  that  treats  heat  conduction,  electromagnetic  field  propagation,  and  temperature-dependent 
absorption.  We  consider  two  specific  cases:  the  temperature-induced  shift  of  the  absorption  peak  of  bound  I2 
excitons  in  CdS  single  crystals  and  the  temperature-induced  shrinkage  of  the  band  gap  in  GaAs/GaAlAs  quantum- 
well  material.  For  both  cases  it  is  shown  that  the  inhomogeneous  spatial  distribution  of  the  nonlinear  absorption 
profile  along  the  direction  of  propagation  of  the  electric  field  undergoes  a  first  -order-like  phase  transition  commen¬ 
surate  with  the  transition  of  the  output  intensity  from  one  value  to  another  for  fixed  input  intensity.  Stability 
conditions  are  presented  for  the  steady  states  of  the  output  intensity  as  a  function  of  the  input  intensity  for  the  cases 
treated.  There  is  a  multiplicity  of  bistability  curves,  and  the  stability  analysis  produced  the  expected  branches  of 
stable  solutions. 


$1 


1.  INTRODUCTION 

There  has  been  a  recent  surge  of  interest  in  intrinsic  optical 
bistability1  (IOB),  i.e.,  optical  bistability  that  does  not  de¬ 
pend  on  external  feedback  such  as  from  the  mirrors  of  a 
high-Q  cavity.  Since  the  first  predictions  for  the  conditions 
for  IOB,2:l  other  theoretical  calculations  have  expanded  on 
the  phenomenon  for  a  wide  variety  of  circumstances.4-17 
The  first  reported  experimental  study  of  IOB  was  conducted 
by  Hajto  and  Janossy18  in  amorphous  GeSe2  due  to  thermal¬ 
ly  induced  absorption,  and  this  was  followed  by  the  experi¬ 
mental  observations  of  Miller  et  al.19  in  GaAs/GaAlAs  by 
tuning  just  below  the  band  gap.  They  also  analyzed  their 
results  in  terms  of  temperature-induced  absorption.  Boh- 
nert  et  a l.,20  and  almost  simultaneously  Rossmann  et  al.,2' 
observed  IOB  in  CdS.  The  presence  of  IOB  in  CdS  has  been 
interpreted  as  due  primarily  to  effects  of  electronic  origin 
owing  to  saturation  of  absorption  because  of  the  generation 
of  carriers  resulting  in  nonlinear  renormalization  of  the  band 
gap  with  increasing  electromagnetic  field  intensity. 1(1  Later, 
Dagenais  and  Sharfin  observed  IOB  in  CdS  at  much  lower 
intensities  by  tuning  just  below  the  bound  I2  exciton  reso¬ 
nance”  and  interpreted  their  results  as  due  t  Hu  Anally 
induced  increasing  absorption.  Since  these  initial  experi¬ 
mental  studies  were  reported,  there  have  been  many  other 
reported  results  of  observations  of  IOB  in  various  materials 


with  different  mechanisms28  24  as  well  as  predictions  of  IOB 
under  other  novel  conditions.  For  instance,  IOB  with  atoms 
densely  packed  so  that  dipole-dipole  interactions  become 
important  has  been  analyzed  by  Bowden  and  co-work- 
ersii.i  >.is-i7  not  y e t  experimentally  studied. 

In  this  paper  we  analyze  nonlinear  longitudinal  spatial 
inhomogeneit  ies  associated  with  IOB  due  to  induced  absorp¬ 
tion.  We  treat  two  specific  cases:  induced  absorption, 
which  results  when  the  incident  laser  field  is  tuned  just 
below  the  bound-exciton  resonance  of  the  I2  bound  exciton 
in  thin,  uncoated  platelets  of  CdS  single  crystals,  and  in¬ 
duced  absorption,  which  results  from  tuning  the  laser  fre¬ 
quency  below  the  band  edge  in  GaAs/GaAlAs  quantum-well 
structures.  For  both  cases,  the  plane-wave  propagation  of 
the  electromagnetic  field,  coupled  with  heat  conduction  and 
the  associated  boundary  conditions,  is  used  to  derive  the 
steady-state  longitudinal  electric  field  amplitude,  intensity, 
and  temperature  and  the  temperature-dependent  absorp¬ 
tion  distributions  in  the  nonlinear  material.  It  is  shown  for 
each  case  that  the  longitudinal  spatial  absorption  distribu¬ 
tion  in  the  material  undergoes  a  first-order  spatial  phase 
transition  commensurate  with  the  transition  of  the  system 
from  one  bistable  output  intensity  to  another  for  fixed  inci¬ 
dent  intensity.  The  longitudinal  spatial  first-order  phase 
transition  is  generic  to  IOB  and  does  not  exist  for  optical 
bistability  (OB)  in  the  presence  of  an  optical  cavity  or  an 
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etalon  of  high  Q.  In  fact,  it  rounds  out  the  concept  of  IOB  as 
an  example  of  a  first-order  phase  transition  far  from  thermo¬ 
dynamic  equilibrium.  The  first  analysis  of  the  longitudinal 
spatial  phase  transition  for  a  dense  two-level  atomic  system 
was  published  by  Ben-Aryeh  and  co-workers.15-17  With  the 
exception  of  the  recent  work  by  the  present  authors,  all 
previous  treatments  have  ignored  the  spatial  distribution 
dependence  of  the  field  and  temperature  in  IOB  due  to 
induced  absorption. 

It  should  be  noted  here  that  the  predictions  of  Koch  et 
air 5  and  subsequent  experimental  observations  of  Gibbs  et 
a l.26  of  the  stepwise  output  as  a  function  of  pulsed  input 
intensity  for  CdS,  which  they  describe  as  “kinks,”  is  not  at 
all  the  same  phenomenon  that  we  describe  here.  In  that 
case,  the  sample  of  CdS  (or  CdS-doped  glass)  is  excited  by  a 
pulse;  an  instability  in  the  transmission  is  found  to  be  due  to 
band-gap  renormalization.  What  we  discuss  for  IOB  here  is 
a  genuine,  stationary,  first-order  spatial  phase  transition. 

In  addition  to  the  treatment  of  the  longitudinal  spatial 
phase  transition,  we  discuss  the  stability  conditions  of  the 
stationary  states  for  the  cases  presented.  We  develop  here, 
for  the  first  time  to  our  knowledge,  a  linearized  stability 
analysis  that  takes  into  account  the  spatial  dependence  and 
intrinsic  longitudinal  inhomogeneity  in  IOB. 

The  model  is  presented  and  discussed  in  Section  2  in 
terms  of  coupled  constitutive  Maxwell  and  heat-conductiv¬ 
ity  equations.  Section  3  is  used  to  present  the  results  of 
numerical  calculations  for  IOB  in  CdS  and  GaAs/GaAlAs 
quantum-well  material.  The  linearized  stability  analysis  is 
developed  in  Section  4,  which  takes  into  account  propaga¬ 
tion  and  spatial  inhomogeneity.  This  analysis  is  applied  to 
the  cases  treated,  and  the  results  are  discussed.  The  final 
section  is  used  for  conclusions  and  to  summarize  the  results. 


2.  THE  MODEL 

Here  we  review  the  model  presented  earlier.1114  Our  model 
for  IOB  due  to  induced  absorption  consists  of  a  set  of  cou¬ 
pled  Maxwell  equations  for  the  electric  field  amplitude  in¬ 
side  the  nonlinear  medium  and  heat-conductivity  equation 
for  the  temperature  distribution  in  the  material.  In  the 
plane-wave  limit,  the  positive-frequency  part  of  the  field, 
7?+,  in  the  medium  is  given  by 

E\x,  t)  =  [E+F(x,  t)eikx  -I-  E+r(x,  t)e~ikx]e~'“* ,  (1) 

where  E(x,  t)  =  E*(x,  t)  +  E~(x,  t ),  (E+)*  =  E~,  and  Ep  and 
Ep  are  forward-  and  backward-propagating  electric  field 
amplitudes,  respectively,  and  are  assumed  to  be  slowly  vary¬ 
ing  functions  of  x  and  t.  Here  k  =  knJ( ,  where  k0  is  the  wave 
vector  in  vacuum  and  t  is  the  real  part  of  the  dielectric 
function  in  the  medium. 

For  steady  state  and  in  the  limit  of  fast  relaxation  of  the 
medium.  Maxwell’s  equations  in  the  slnwlv  varying  envelope 
approximation  (SVEA)  are 

1  dEp  dEp  „  x 

i^f  +  ^  +  l£P  =  0’  (2a) 

c  dt  dx  2 

1  dEt  dEt  „ 

--'ht  +  it-oE+b- o,  (2b) 

c  dt  dx  2 

where  n  =  a  +  10  and  a  and  0  are  corresponding  absorption 
and  dispersion  functions  in  the  medium  and  are  related  by 
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the  Kramers-Kronig  relations.  It  is  to  be  noted  that  r)  is 
spatially  inhomogeneous  in  the  medium  and  is  dependen*  on 
the  local  temperature  T.  Thus  Eqs.  (2)  are  coupled  by  rj. 
The  temperature  Tlx)  is  determined  by  the  heat-conduction 
equation: 

(3) 

‘  dt  dxt 

where  k  is  the  heat  conductivity,  assumed  constant  and  ho¬ 
mogeneous,  and  I  =  l(x)  is  the  local  intensity  of  the  field.  In 
this  section  and  in  Section  3  we  consider  only  steady-state 
solutions,  but  the  time  dependence  is  important  in  Section  4, 
where  we  take  up  the  stability  analysis. 

Maxwell’s  equations,  Eqs.  (2),  can  be  simplified  further  by 
using  the  transformation 

Epr,(x)  =  E^UJexpj^  i  |  dx'd(x')  •  (4) 

The  resulting  equations  thus  have  real  coefficients,  and  the 
dispersion  0  appears  only  in  the  boundary  conditions  for  the 
field  at  x  =  L  in  terms  of  a  phase  factor  denoted  by 

—  If L 

KL  =  KL  -  —  dx'0(x').  (5) 

2  Jo 

Consistent  with  the  SVEA,  the  intensity  I  in  Eq.  (3)  can  be 
written  in  terms  of  the  forward-  and  backward-propagating 
fields  as 

/  =  I  Ep\2  +  \Ep\2,  (6) 

where  the  rapidly  varying  cross  term  has  been  neglected. 
Equations  (2),  using  the  transformation  [Eq.  (4)],  become 

dEp  a  _x 

ir+2^  °’  (7a) 

dEfi  CY  — 4. 

aT-2^"°-  (7b) 

and  (refer  to  Fig.  1)  the  associated  boundary  conditions  are 


at  r  =  0 


and  transmitted  field  amplitudes,  respectively,  and  (()  is  the 
dielectric  constant  of  the  external  medium,  whereas  e  is  the 
linear  dielectric  constant  of  the  medium. 

In  terms  of  the  incident  intensity  7/  and  the  transmitted 
intensity  It,  the  boundary  conditions,  Eqs.  (8),  give 


h  [l  +  (rA0)1/2l2ML)’  (9> 

I,  =  \a2IF(0)  +  b2IH(0)  +  2a b  Re[Ef+(0)£fl(0)]|,  (10) 


JHC 

■Vi 

m 


w 

m 

m 


wv 
1*1, ,‘« 

Mil 


Ep(  0)  +  E*h{  0)  =  EJ  +  Ep, 

(8a) 

f&Si 

v(l 

i<[EpW)  ~  E+H( 0)|  =  (EJ  -  Ej)v/q„ 

(8b) 

-  •  J 

■v-y. 

at  x  =  L 

•y-y 

Ep(L)e‘KI-  +  Ep(L)e~'K'-  =  E+r, 

(8c) 

■ 

v  y 

MEtiL)e'*'-  -  Et(L)e-Rl  ]  =  E+Tfu, 

(8d) 

where  £/,  7?r,  and  Ej  are  the  normally  incident,  reflected. 

i 

H 


Bowden  el  at. 


Vol.  5,  No.  l/.lanuarv  198H/.I.  Opt.  Sue.  Am.  B 


where 

o  =  (f/t,,)1'2  +  1,  b  -  -  1,  (11) 

and 

l  fix)  =  \Et(x)\-,  IH(x)  =  \E+h{x)\'\ 

I(x)  =  \Et-(x)\2 +  \EJi{x)\--  (12) 

For  the  boundary  conditions  associated  with  Eq.  (3),  we 
have  adopted  the  following  form: 

Hx\T{x)-Th\  =  HxMx)XxZ°r  ( 13) 

where  and  H /  are  related  to  the  geometry  and  the  thermal 
resistance  at  the  respective  boundaries  and  7*  is  the  ambient 
temperature.  Since  a  =  nfTU)]  is  a  nonlinear  function  of 
the  temperature  in  the  medium,  Eq.  (13)  constitutes  a  non¬ 
linear  relationship  between  the  temperature  T  at  each 
boundary  and  the  field  intensity  /  at  the  respective  bound¬ 
ary.  Thus  the  local  bistable  condition  throughout  the  medi¬ 
um  is  controlled  at  the  boundaries. 

The  set  of  coupled  constitutive  Maxwell  equations  and  the 
heat-conductivity  equation,  Eqs.  (2)  and  (3),  together  with 
their  associated  boundary  conditions,  Eqs.  (8)  and  (13),  rep¬ 
resent  our  model  for  IOB  due  to  induced  absorption.  Here 
we  apply  this  model  to  two  distinct  cases,  which  we  treated 
previously,  to  discuss  a  new  aspect  of  the  first-order  phase 
transition,  namely,  the  longitudinal  spatial  inhomogeneity 
associated  with  the  intrinsic  local  bistability  in  the  medium. 
The  two  cases  that  we  shall  address  are  the  following:  Case 
A.  the  situation  in  which  the  incident  field  is  tuned,  at  low 
field  values,  just  below  the  absorption  resonance  for  the 
bound  f.  exciton  in  CdS  single-crystal  platelets;  and  Case  B, 
in  which  the  laser  field  is  tuned  just  below  the  band  edge  in  a 
GaAs/GaAlAs  quantum-well  structure.  With  regard  to 
Case  A,  the  change  in  the  nonlinear  absorption  as  a  function 
of  the  internal  field  intensity  causes  nonlinear  inhomogen¬ 
eous  heating  of  the  material  with  a  consequent  local  shift  of 
the  absorption  peak  to  lower  values  with  increasing  local 
temperature;  whereas  in  Case  B,  internal  heating  of  the 
material  causes  shrinkage  of  the  band  gap  to  lower  values  for 
increasing  local  temperature. 

The  model  presented  here  was  used  previously  to  analyze 
both  cases.  The  absorption  «  and  the  dispersion  d,  as  func¬ 
tions  of  local  temperature  T(x),  are,  however,  of  entirely 
different  origin  in  the  two  cases  and  are  of  completely  differ¬ 
ent  functional  form.  The  boundary  conditions  are  also 
quite  different  in  the  two  cases.  Each  case,  however,  repre¬ 
sents  a  generic  prototype,  as  will  become  apparent  in  the 
development  in  the  next  two  sections.  In  Section  3  we 
discuss  results  for  longitudinal  spatial  inhomogeneities  and 
first-order  phase  transitions  corresponding  to  IOB  due  to 
induced  absorption  in  the  two  kinds  of  system. 

3.  LONGITUDINAL  SPATIAL  FIRST-ORDER 
PHASE  TRANSITIONS 

A.  Case  A:  CdS 

Here  we  discuss  IOB  due  to  induced  absorption  associated 
with  the  resonant  absorption  nonlinearity  in  CdS  due  to  the 
bound  Ij  exciton.  As  previously  discussed,1"  we  adopt  here 
the  empirical  form  for  the  function  ij  of  Eqs.  (2): 


i:t 


n(x) 


o„|7'(-r)) 

1  +  ok  -  w0[7’(x)|| 


(14) 


with  the  explicit  form  for  the  absorption  part  of  Eq.  (141 
taken  from  experimental  data  for  77,  =  7  K: 


1  |3  +  0.04(7'-  r,,)-] 


10  |l  + 


„m|-(().4>-| 


(15) 


in  units  of  reciprocal  micrometers,  where  T  is  in  degrees 
Kelvin  and  (w  —  w„)  is  in  the  unit  of  reciprocal  centimeters. 
In  addition,  the  following  empirical  relation  is  used  for  the 
temperature,  T,  dependence  of  the  peak  of  the  absoiptioii. 


a !„(T)  =  a,'„(0)  -  AT1. 


(16) 


This  relation  is  valid  for  0  <  T  <  40  K  with  the  parameter 
values  A  =  1/80  cm'1  K'-  and  u.\i(0)  =  20531  cm'1.  The 
conductivity  *  of  Eq.  (3)  is  taken  as  k  =  I0'4  W /pm  K. 

As  is  discussed  in  Ref.  13,  the  IOB  in  this  case  is  extremely 
sensitive  to  the  heat-conductivity  boundary  conditions,  Eq. 
(13)  and,  in  fact,  is  dependent  or.  these  conditions,  whereas 
the  boundary  conditions,  Eqs.  (8)  for  the  field  equation^, 
Eqs.  (7),  do  not  cause  bistability.  We  discussed  previously 
that  the  threshold  conditions  and  hysteresis  in  the  output 
intensity  as  a  function  of  the  input  intensity  /„,  can  be 
strongly  affected  by  changing  the  boundary  conditions  asso¬ 
ciated  with  the  heat  conduction. 

In  a  manner  similar  to  that  discussed  in  Ref.  13,  we  have 
numerically  integrated  Eqs.  (7)  and  (3),  using  the  associated 
boundary  conditions,  Eqs.  (8)  and  (13),  and  the  empirical 
relations  (15)  and  (16)  for  T>,  =  0  and,  for  Eq.  (13),  Hi.  =  HJ 
3.  Because  the  heat-conductivity  boundary  conditions  are 
suitably  chosen  at  the  two  ends  of  the  CdS  platelet  (Fig.  1), 
we  obtain  two  hysteresis  conditions,  one  due  to  increasing 
induced  absorption  at  lower-field  values  and  the  other  due 
to  decreasing  induced  absorption  at  the  higher  fields.  The 
results  are  qualitatively  the  same  as  those  reported  earlier 
and  are  shown  in  Fig.  2.  The  only  difference  between  the 
results  shown  in  Fig.  2  and  those  reported  earlier1'  is  that  in 
this  case  we  have  chosen  the  ambient  temperature  77,  =  0  for 
convenience. 

The  absorption  profile  <v(jf)  before  switching  from  the 
high  to  the  low  transmission  state  is  shown  in  curve  (a)  of 
Fig.  3,  which  corresponds  to  point  a  in  the  inset  of  Fig.  2. 
After  switching  to  the  low  transmission  state,  for  the  same 
input-field  intensity,  the  absorption  profile  «(.v)  changes  to 
that  shown  in  (b)  of  Fig.  3,  which  corresponds  to  point  b  of 


T(O) 


T(L) 


Fig.  1.  CdS  configuration. 
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Fig.  2.  CdS  increasing  and  decreasing  absorption  hysteresis  in 
scaled  intensities,  J  =  I/H:  b,  lower  field  (increasing  absorption 
regime);  c,  higher  field  (decreasing  absorption  regime). 


Fig.  3.  «(tr)  for  the  case  of  Fig.  2. 

Fig.  2.  Thus  the  peak  of  the  absorption  has  swept  through 
the  fixed  field  frequency  in  the  medium  for  a  portion  of  the 
nonlinear  medium  nearest  the  left-hand  boundary,  and,  con¬ 
sequently,  the  peak  in  the  spatial  distribution  of  the  absorp¬ 
tion  appears  in  the  medium,  just  inside  the  input  boundary. 
The  upper  and  lower  curves  for  «(x)  in  Fig.  3  are  the  spatial 
absorption  distribution  profiles  corresponding  to  the  low 
and  high  transmission  states  for  the  same  value  of  the 
threshold  input  field  intensity  [Fig.  2  (inset)].  Thus  the 
spatial  distribution  for  the  absorption  «(x)  changes  abrupt¬ 
ly,  corresponding  to  the  switching  of  the  output  intensity 
from  high  to  low  transmission,  and  constitutes  a  manifesta¬ 
tion  of  the  first-order  phase  transition.  This  aspect  of  the 
phase  transition  is  intrinsic  to  IOB  and,  of  course,  cannot 
occur  in  optical  bistability  that  depends  on  optical  feedback, 
as  in  a  nonlinear  Fabry-Perot  or  ring  cavity. 

The  corresponding  intensity  distributions  in  the  medium, 
/(x),  are  depicted  in  Fig.  4,  and  the  corresponding  slowly 
varying  forward-  and  backward-propagating  field  ampli¬ 
tudes  are  shown  in  Fig.  5.  Also  noted  is  the  change  in  the 
overall,  or  average,  temperature  gradient  across  the  sample 
by  more  than  a  factor  of  10  associated  with  the  transition. 
In  the  change  to  the  lower  transmission  state,  the  tempera¬ 
ture  in  the  sample  at  the  input  boundary  has  increased  by 
13.5  K  because  of  increased  absorption,  whereas  the  corre¬ 
sponding  temperature  at  the  output  boundary  has  dropped 
by  3.1  K,  which  is  due  to  the  drastic  decrease  in  field  intensi¬ 
ty,  7(x),  near  that  boundary. 

The  heat-conductivity  boundary  condition,  Eq.  (13),  us¬ 
ing  Eqs.  (15)  and  (16),  is  shown  in  Fig.  6  for  the  input  side. 
This  shows  that  the  low-field  first-order  phase  transition 


discussed  above  is  primarily  associated  with  the  input-side 
thermal  boundary  condition,  and  the  qualitative  aspect  of 
the  behavior  depicted  in  Figs.  2-5  is  dependent  on  this 
boundary  condition.  Since  there  is  a  boundary  condition  of 
the  same  nature  associated  with  the  output  boundary,  and 
since,  as  air  ady  noted,  the  temperature  at  that  boundary 
has  decreased  in  the  transition  from  the  high  to  low  trans¬ 
mission  states,  it  is  anticipated  that  another  phase  transition 
will  occur  at  higher  input  intensity.  This  transition  is 
shown  in  Fig.  2  and  is  due  to  decreasing  absorption. 

Thus,  as  the  incident  intensity  is  increased,  the  peak  of  the 
absorption  profile,  o(x),  appears  at  larger  values  of  x  owing 
to  the  monotonic  increase  in  the  intensity  in  the  medium, 
/(x),  and  the  distribution  «(x)  is  reshaped  longitudinally. 
The  behavior  for  o(x)  is  displayed  in  Fig.  7  for  two  different 
input-field  values  as  indicated  in  Fig.  2  along  the  lower 
output-field  branch.  In  other  words,  as  the  internal  intensi¬ 
ty  /(x)  is  increased  further  throughout  the  medium,  the  local 
absorption  peak,  wolTXx)],  shifts  to  lower  values,  and  since 
the  laser  frequency  oij.  was  set  at  a  value  less  than  o>(>  at  low 
input  fields,  the  absorption  resonance  is  swept  through  the 
laser  frequency  for  a  larger  portion  of  the  material  medium. 
This  behavior  is  evident  from  Eqs.  (15)  and  (16). 

The  local  IOB  of  the  material  is  expressed  explicitly  in 
terms  of  the  boundary  condition  (Fig.  6).  However,  because 
of  the  large  and  dominant  diffusion  in  this  case,  because  the 
heat  conductivity  k  is  large  [Eq.  (3)],  the  spatial  phase  transi¬ 
tion  corresponds  to  a  kind  of  Maxwell  construction  far  from 
thermodynamic  equilibrium.  Thus  the  intensity  distribu- 


Fig.  4.  /(x)  for  conditions  corresponding  to  Fig.  3:  (a)  before 

switching;  (b)  after  switching. 


Fig.  5.  (a)  £/j(x)  and  £f(x)  before  switching  (Fig.  3).  (b)£«(x)and 

Ey(x)  after  switching  (Fig.  3). 
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Fig.  6.  AT  versus  J\„  at  the  input  face  from  Eq.  (13)  for  conditions 
of  Fig.  3.  Here,  Jjn  =  I/H. 


Fig.  7.  Input  field  /independence  of  the  nonlinear  absorption,  «(*). 

tion  [/(*),  Fig.  4]  is  monotonic.  For  vanishingly  small  diffu¬ 
sion  in  a  transport  equation,  we  would  expect  an  abrupt 
change  in  Hx),  at  some  0  <  x  <  L,  introducing  a  spatial 
discontinuity  in  the  nonlinear  index  An  and  a  sharp  bound¬ 
ary  in  the  material  between  two  spatially  separated  phases. 
This  kind  of  phase  transition  was  referred  to  earlier.17  The 
case  treated  here  is  for  the  opposite  limit,  i.e.,  large  diffusion, 
which  is  the  natural  situation  for  thermally  induced  absorp¬ 
tion. 

A  second,  or  high-field,  threshold  is  reached,  which  is 
shown  in  Fig.  2  for  the  output  intensity  versus  the  input 
intensity  Im.  The  spatial  absorption  profiles  at  threshold, 
corresponding  to  the  low  and  high  transmission  states  for 
the  same  threshold  input  field,  are  shown  in  Fig.  8.  In  this 
case,  the  absorption  peak  is  switched  out  of  the  medium  as 
the  input  field  is  increased,  and,  in  contrast  to  the  low-field 
phase  transition,  this  situation  corresponds  to  switching  due 
to  decreasing  absorption.  Thus  the  hysteresis  loop  is  coun¬ 
terclockwise,  as  opposed  to  the  clockwise  hysteresis  associat¬ 
ed  with  the  low-input-field  transition,  which  is  due  to  in¬ 
creasing  absorption. 


It  is  noted  that  the  temperature  at  the  output  boundary 
has  increased  by  nearly  a  factor  of  4,  where,  in  contrast,  the 
temperature  at  the  input  boundary  remains  nearly  fixed 
with  respect  to  the  transition.  Thus  the  transition  dis¬ 
cussed  here  is  associated  primarily  with  the  output  bound¬ 
ary  condition,  Eq.  (13).  The  large  increase  in  the  tempera¬ 
ture  at  the  ouLput  boundary  is  due  to  the  relatively  large 
increase  in  the  intensity  l(x)  near  the  output  boundary. 
Here,  the  absolute  value  for  the  average  temperature  gradi¬ 
ent  across  the  sample  drops  in  the  transition  from  the  low  to 
the  high  transmission  state.  Noted  here  is  an  actual  inver¬ 
sion  of  the  gradient,  so  the  temperature  at  the  output  is 
greater  than  the  temperature  at  the  input,  i.e.,  T(L)  >  T( 0). 

If  the  input  field  is  decreased  from  a  h'gh  value,  an  inter¬ 
mediate  field  threshold  is  reached,  as  depicted  in  Fig  2. 
Figure  9  shows  the  absorption  profile  a(x)  for  the  high  and 
low  transmission  states,  respectively,  at  the  indicated 
threshold  input  field.  Here  the  peak  of  the  absorption  pro¬ 
file  o(x)  switches  back  into  the  material  in  the  transition 
from  the  high  to  the  low  transmission  state  (Fig.  9).  It  is 
interesting  to  compare  the  profile,  o(x),  in  the  low  transmis¬ 
sion  states  in  Figs.  9  and  8.  In  the  former  case,  the  absorp¬ 
tion  peak  is  located  deeper  into  the  material  as  compared 
with  the  latter  case,  which  is,  of  course,  expected,  owing  to 
the  difference  in  the  respective  input  fields. 

We  note,  in  connection  with  Fig.  9,  that  T(0)  remains 
approximately  fixed  but  T(L)  decreases  markedly  through 
the  transition,  indicating  that  the  phase  transition  associat¬ 
ed  with  the  high-field,  counterclockwise  hysteresis  is  associ- 


Fig.  8.  «(x)  at  higher  input  field,  second  threshold:  (e)  low  trans¬ 
mission  (Fig.  2);  (I)  high  transmission  (Fig.  2). 


Fig.  9.  i*(x)  at  lower  input-field  threshold  in  high-input-field  re 

gime:  (g)  high  transmission  (Fig.  2);  (h)  low  transmission  (Fig.  2). 
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Fig.  10.  o(.v)  at  lowest  input-field  threshold:  (i)  low  transmission 

(Fig.  2);  (i)  high  transmission  (Fig.  2). 


ated  primarily  with  the  output-side  heat -conductivity 
boundary  conditions. 

Finally,  for  still  lower  values  of  the  input-field  intensity, 
the  lowest  field  threshold  is  reached,  as  shown  in  Fig.  2. 
Figure  10  shows  the  absorption  profiles  o(.r)  in  the  low  and 
high  transmission  states,  respectively.  Here,  the  tempera¬ 
ture  T(L)  changes  relatively  little,  indicating  that  the  transi¬ 
tion  is  controlled  primarily  by  the  heat-conductivity  bound¬ 
ary  conditions  at  the  input  boundary. 

B.  Case  B:  GaAs/GaAlAs 

For  Case  B  we  discuss  10B  due  to  induced  absorption  using 
the  nonlinear  absorption  associated  with  tuning  in  the  re¬ 
gion  of  (he  band  gap  in  GaAs/GaAJAs  multiple-quantum- 
well  structures.  Here  we  use  the  empirically  determined 
absorption  function 


T -(w  -  »,)-]  C -(«  -  u»v-*i 

u(u>,  7’)  =  <»,  exp  2r,  I  +  »/,  expl  i)(^  i 

+  2,.,[l+exp(“^W)J 

X  |1  +  exp<— 27r/?v,'2|S2,  -  wlM/-)r',  (17) 

where  all  the  constants,  <«/,  F/,  <>/,,  17, ,  II/,  U/„  II  ,  F, ,  and  Rv, 
are  given  in  Ref.  14.  The  band  gap  has  the  temperature,  77 
dependence. 

AT - 

12,(77  =  12,(0)  -  ■  (18) 

/  +  n 

where  !!.(<))  =  1.51216  eV,  A  =  8.87  X  lO"4  K1,  H  =  572  K, 
and  T  is  in  units  of  K.  The  results  discussed  here  pertain 
exactly  and  explicitly  to  the  calculations  and  results  report¬ 
ed  earlier  and  appear  in  Ref.  14.  Equation  (17),  o(77  versus 
77  is  shown  in  Fig.  1 1 . 

The  geometry  of  the  configuration  is  discussed  in  Ref.  14. 
Equations  (2)  and  (3)  are  integrated  across  the  sample  under 
steady-state  conditions  and  with  the  appropriate  boundary 
conditions,  as  discussed  in  Ref.  14.  The  heat-conductivity 
boundary  conditions,  Eq.  (13),  are  shown  for  this  case  in 
Figs.  12(a)  and  12(b)  for  the  input  boundary  and  the  output 
boundary  of  the  sample,  respectively.  The  corresponding 
output  intensity  versus  the  input  intensity  /,„  is  repre¬ 
sented  in  Fig.  13.  These  results  are  identical  to  the  results 
reported  earlier  but  are  reproduced  here  for  clarity. 


The  absorption  profile,  o(x),  is  shown  in  Fig.  14  for  six  of 
the  eight  threshold  conditions  depicted  in  Figs.  12  and  13. 
At  the  low-input -field  threshold,  the  transition  from  the 
high  to  the  low  transmission  state,  denoted  by  a  and  I), 
respectively,  in  Figs.  12  and  13,  corresponds  to  the  change  in 
the  absorption  distribution,  oU).  in  the  material  denoted  by 
curves  a  and  b  in  Fig.  14.  Thus  the  absorption  profile  of  the 
band  edge  switches  into  the  medium  as  the  sample  switches 
from  the  high  to  the  low  transmission  state.  The  corre¬ 
sponding  transition  at  the  input  boundary  is  noted  in  Fig.  12, 
and  it  is  seen  that  the  transition  is  associated  primarily  with 


2B0  300  3J0  340  360  380  400 


Fig.  1 1 .  n(  7’)  versus  T  for  (iaAs/GaAIAs. 


Fig.  1 2.  (a)  77,i  versus  /,„  for  GaAs/GaAIAs.  (1>)  T„M  versus  for 
GaAs/GuAIAs. 


Fig.  13.  versus /,„  for  conditions  of  Figs.  11  and  12. 
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Fig.  14.  aix)  for  GaAs/GaAlAs.  The  conditions  are  the  same  as  for 
Figs.  11-13. 


the  input  boundary,  whereas  the  corresponding  changes  at 
the  output  boundary  are  relatively  minor. 

For  higher  values  of  the  incident  field  intensity,  the  edge 
of  the  absorption  profile  penetrates  deeper  into  the  medium 
until  the  higher  field  threshold  is  reached,  denoted  by  e-f  of 
Figs.  12  and  13.  The  corresponding  absorption  profiles  a(x) 
for  the  high-  and  low-transmission  states  for  this  threshold 
condition  are  shown  in  Fig.  14  and  are  denoted  by  e  and  f, 
respectively.  Thus  the  edge  of  the  absorption  profile  «(x)  is 
swept  out  of  the  medium.  This  transition  is  controlled  pri¬ 
marily  by  the  boundary  conditions  at  the  output  boundary, 
as  noted  in  Fig.  12. 

Unlike  the  situation  for  Case  A,  for  which  «(x)  has  a  well- 
defined  absorption  maximum,  the  high-  and  low-field  tran¬ 
sitions  each  correspond  to  transitions  due  to  increasing  ab¬ 
sorption  and  thus  give  rise  to  clockwise  hysteresis  loops.  In 
other  words,  as  the  incident  field  intensity  is  increased,  the 
absorption  distribution  «(x)  in  the  material  increases  mono- 
tonically  to  saturation. 

In  this  section  we  have  treated  the  longitudinal  spatial 
first-order  phase  transition  associated  with  the  spatial  de¬ 
pendence  of  the  absorption  distribution  «(x)  in  the  material 
for  IOB  due  to  induced  absorption.  This  gives  rise  to  a  kind 
of  Maxwell  construction  in  the  material  far  from  thermody¬ 
namic  equilibrium  and  rounds  out  the  profile  of  a  first-order 
phase  transition  interpretation  of  IOB.  This  bears  out  the 
intrinsic  nature  of  the  phase  transition  in  IOB. 

In  Section  4  we  discuss  the  stability  conditions  for  the 
cases  treated  in  this  section. 

4.  STABILITY  ANALYSIS 

We  shall  address  the  stability  analysis  to  the  CdS  system; 
the  identical  procedure  is  applicable  to  the  GaAs/GaAlAs 
system  treated  here,  and  the  results  are  qualitatively  the 
same.  Considering  the  general  solution  of  the  physical 
quantities  such  as  Epjt(x,  t)  as  the  sum  of  the  steady-state 
solution  (£f  fl(x)],  obtained  earlier  and  a  small  transient 
variation,  bEppix,  t ),  i.e., 

Efh(x,  t )  =  [£f  W(x)]s  +  6EFH(x,  t ), 

T(x,t)  =  [r(x,t)K  +  &T(x,t), 

we  obtain  from  Eqs.  (2) 


+  Je,,,  lar=  o  09) 


C,  8&J  ~  *  a'„  «T  =  /f“  ,  AT  +  2a(Ef)JtEF  +  2«<E,., ).«£„, 

Of  dx~  al 


where  (  >,  denotes  the  steady  state  and  we  have  dropped  the 
superscript  +  for  convenience.  The  usual  stability  analysis 
is  carried  out  by  solving  the  linearized  differential  equations 
and  fetermining  a  set  of  eigenvalues  j\!  that  express  the  time 
development  through  the  relations 

SEFM(x,t)  =  b£pH{x)eX1,  bT(x,t)  =  bT(x)e".  (21) 

The  numerical  solution  corresponding  to  this  procedure  has 
been  carried  out  by  one  of  the  authors-'  using  finite-differ¬ 
ence  methods.  We  present  here  a  different  approach  based 
on  approximations  intrinsic  to  the  particular  systems  treat¬ 
ed  here.  This  approach,  though  less  general  than  the  con¬ 
ventional  one,  has  the  advantage  of  producing  transparent 
analytical  results.  The  numerical  results27  are  in  agreement 
with  those  given  here. 

Because  of  the  high  diffusion  limit  applied  to  the  cases 
treated  here,  the  heat-conduction  equation,  Eq.  (3),  cannot 
contribute  to  instabilities  associated  with  the  steady  state. 
This  is  evidenced  by  the  fact  that  the  temperature  gradient 
in  the  material  is  essentially  linear.  Furthermore,  as  we 
have  demonstrated  in  the  previous  sections,  the  bistable 
behavior  of  the  system  is  entirely  associated  with  the  multi¬ 
valued  solutions  of  the  boundary  values  at  the  two  ends  of 
the  sample  [Eq.  (13)].  Consequently,  under  the  approxima¬ 
tion  that  bistable  switching  is  independently  controlled  by 
the  nonlinear  boundary  conditions  at  the  two  end  faces  of 
the  material,  we  analyze  the  stability  at  ,r  =  0  and  at  x  =  L 
separately. 

First,  consider  Eq.  (19)  evaluated  at  x  =  0: 

-  bEFH(x)  ±  ~  bEp /f[x)  +  ”  ,  bEpi,(x) 
c  dx  2 


AT(.r)  =  0;  (22) 


and  AT( 0)  and  bEpn(x)  at  x  =  0  are  related  through  Eqs.  (8) 
and  (13): 


.  (9cy  i 

bT(x)  =  a  1  -  -r-rI(x)  [E„ (x)  +  REF{ j 
O I  J  s 


x )] «£„(*). 


bEH(x)  =  bEp(x),  (23b) 

n 

where  the  intensity  I  is  normalized  by  taking  the  surface 
conductivity  to  be  unity  and  R  =  (Vf  -  l)/(yr  +  1). 

If  we  express  bEp  n(x)  in  terms  of  their  Fourier  transforms, 

bEp H(x)  =  'S  exp(±2jr inx/L)bEp (24) 


18  J.  Opt.  Soc.  Am.  B/Vol.  5,  No.  1/January  1988 


Bowden  ft  al. 


where  the  upper  (lower)  sign  is  taken  for  HEy  (bEp),  and 
eliminate  bT(x),  using  Eq.  (23a),  in  Eqs.  (22),  these  equa 
tions  yield 


(A  +  iknc  +  gr)bEF["'  -  0, 
(A  -  iknc  +  gH)bEHln'  =  0, 


where,  at  x  =  0, 


gF  ~  «(0) 


gH  =  a(0)  — 


+  2(~)  eb(0)eju<  - 1) 

\dJ  /x=0 

j  _  /d<v\  -  rjr 

l  \<9T  }x=o  a(0)  J 

1-2(ILo^0)£-"/(^-1) 

r _ / do \  Trm ->»]}' 

j  \dT  )x=o  a(0)  J 


(25a) 

(25b) 


(26a) 


(26b) 


Notice  that  if  Eq.  (24)  and  a  similar  expansion  for  bT(x)  are 
used  in  Eqs.  (19)  and  (20),  the  conventional  approach  to 
instabilities  emerges,  whereas  replacing  the  physical  quanti¬ 
ties  nl,„  oEflL  etc.  by  average  values  in  the  medium  leads  to 
results  in  the  mean-field  approximation. 

As  is  evident  from  Eqs.  (25)  and  (26),  stability  of  the  phase 
transition  associated  with  the  left-hand  boundary,  i.e.,  x  =  0, 
is  determined  by  the  sign  of  gy  and  gB,  i-e.,  both  gy  and  gn 
must  be  positive  to  warrant  a  stable  solution.  An  analysis  of 
gy  and  gB  shows  that  the  dominant  factor  in  the  determina¬ 
tion  of  the  sign  is 


A()=  1  -  (dn/d'f '),„,(  1  (0)  -  Th ]/«(0).  1 27 ) 

The  reason  for  this  is  that  the  numerator  of  gy  is  close  to  one, 
since  E«(0)  is  small,  whereas  the  numerator  of  gB  is  always 
slightly  greater  than  the  denominator.  Thus  the  stable  so¬ 
lutions  are  given  by 

A, |  >  0  at  x  =  0.  (28) 


It  is  noted  that  A»  =  0  is  the  exact  location  of  the  turning 
point  (Fig.  6)  at  x  =  0.  Similar  results  pertain  to  the  right- 
hand  boundary,  i.e.,  at  x  =  L,  and  the  phase  transition 
associated  with  that  boundary. 

The  stable  solutions  are  defined  as  the  region  where  both 
Im(A„)  are  negative.  Then  the  unstable  region  in  our  prob¬ 
lem  is  simply  the  middle  branch  in  the  OB  curve,  as  expected. 


5.  CONCLUSIONS 

We  have  demonstrated  and  discussed  the  existence  of  a 
spatial  longitudinal  first-order  phase  transition  generic  to 
IOB  due  to  induced  absorption  and  the  associated  heat- 
conductivity  boundary  conditions.  The  phenomenon  is  in¬ 
deed  intrinsic  to  IOB  and  is  nonexistent  in  the  good-cavity 
limit. 

Here,  we  have  addressed  the  case  of  the  high-diffusion 
limit  characterized  by  a  large  Maxwell  construction  region  in 
the  longitudinal  spatial  variable  characteristic  of  IOB  due  to 
induced  absorption.  This  is  contrasted  with  the  low-diffu¬ 
sion  limit  treated  earlier,17  which  exhibits  an  abrupt  change 
in  the  longitudinal  spatial  dependence  of  the  material  vari¬ 
ables  at  the  phase  transition.  This  behavior  rounds  out  the 
picture  of  an  intrinsic  first-order  phase  transition  in  a  sys¬ 


tem  consisting  of  the  radiation  field  interacting  with  matter 
far  from  thermodynamic  equilibrium. 

The  spatial  first-order  phase  transition  may  be  observed 
directly  by  low-intensity  transverse  field  absorption  in  the 
material.  We  anticipate  that  experiments  will  be  conducted 
in  the  future  to  test  these  fundamental  interpretations. 
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